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Abstract. We study the statistics of column-convex lattice animals resulting from 
the stacking of squares on a single or double staircase. We obtain exact expressions 
for the number of animals with a given length and area, their mean length and 
their mean height. These objects are closely related to Fibonacci numbers. On a 
single staircase, the total number of animals with area k is given by the Fibonacci 
number Fj^. 



1. Introduction 

A lattice animal is a cluster of occupied sites on a lattice. In two dimensions, the area 
of an animal is the number of sites belonging the cluster and its perimeter is defined 
as the set of empty first neighbours of ocupied sites. Alternatively, instead of the site 
clusters, one may consider the corresponding clusters of occupied cells on the dual 
lattice, also called polyominoes. 

The enumeration of lattice animals according to their area and/or perimeter is 
a subject of active research in the statistical physics and combinatorics communities 
(see [1] for a recent review). On the physical side, the main interest lies in the close 
connection between lattice animals and the percolation problem [2,3]. 

In the most general case, the counting problem is quite difficult and only some 
bounds on the asymptotic behaviour are known [4]. This led to the introduction of 
restricted classes of animals (Ferrers graphs, convex and/or directed animals) for which 
some exact results could be obtained, mainly in two dimensions (see for example [5- 
12]). 

In the present work, we consider animals resulting from the stacking of squares 
on a single or double staircase, i.e. column-convex (or vertically convex) animals, for 
which the intersection of a vertical line with the perimeter has at most two connected 
components. These animals are closely related to Fibonacci numbers. In the case 
of a single staircase, the correspondance is particularly simple: the area-generating 
function is equal to the Fibonacci number generating function. 

The paper is organized as follows. In section 2, we study the stacking of squares on 
a single staircase. We calculate the number of animals with a given length and area, 
the total number of animals with a given area, their mean length and mean height 
and study the asymptotic behaviours. The same is done in section 3 for animals on 
a double staircase. The results are discussed in section 4. Some technical details are 
given in the appendix. 
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2. Single staircase 

In this section we consider column-convex lattice animals resulting from the stacking 
of k squares on a single staircase as shown in figure 1. Two neighbouring columns are 
connected when they share at least one edge. 

2.1. Number of animals 

In order to count the number Fj^ i of different animals with area k living on I stairs, 
we introduce the generating function 

CO 

Fiz,t)= J2 (2.1) 

It satisfies the relation 

F{z,t) = zt + t-—[l + F{z,t)] (2.2) 

-L Z 

where the first term corresponds to a single square, the factor z'^ /{I — z) in the second 
term is the generating function of a column with at least two squares, needed for the 
animal to eventually continue its growth on the next stair. Hence we have 

oo 

(2.3) 

= EE = EE 

n=0 p ^-^^ I k=l ^ 



zH'. 



By convention, in sums containing binomial coefficients, the range of summation is 
not explicitly indicated. It is automatically determined by the nonvanishing values of 
the binomial coefficients and, for example, < p < n in equation (2.3). 

The identification of the coefficients of z^t^ in equations (2.1) and (2.3) leads to 

F., - (J:;) . (2.4) 

Using the addition/induction relation for the binomial coefficients [13, p 174], 

/c-/ + l\ fk-l\ fk-l\ 

+ , o (2-5) 



l-l J \l-2, 
ones obtains the recursion relation: 

^k+i,i = ^k,i + ^k-i,l-i ■ (2-6) 
The total number of animals with area k is given by 
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Figure 1. Different stacfcings of five squares on a staircase arranged by number of 
occupied stairs. 



i.e. by the kth Fibonacci number (see [13, p 303]) which satisfies the recursion relation 
= Ff^ + Ff^_^ with the initial values Fq = and F^ = 1, according to (2.6) 
and (2.7). 

This relation with the Fibonacci numbersf can be deduced directly by setting t = 1 
in the first line of (2.3) which leads to 



where F{z) is the generating function of the Fibonacci numbers. 

All the animals built from five squares are shown in figure 1, the different columns 
corresponding to the different values of /. 

2.2. Mean length 

The mean length of animals with k is defined as 



oo 




(2.8) 



k=i 



m 



1=1 ^k,i ^k 




(2.9) 



In order to calculate A^, let us introduce the auxiliary generating function 



oo 



dF{z,t) 
dt 



l-z 



A{z) = Y,A,z 



k 



F{z) + zF'^{z) 



z 



F\z) 



(2.10) 



k=i 



According to [13, p 354] 




(2.11) 



f A similar connection between tfie number of column-convex directed animals with area k and the 
Fibonacci numbers with odd indices i<2fe-i has been noticed a long time ago. See [14]. 
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hence the coefficients of A{z) are 



(2) A;(2F,-Ffc_i) + 3F, 



(2.12) 



where we used the recursion relation for the Fibonacci numbers. The symbol [z'']f{z) 
means the coefficient of z'^ in the series f{z). The mean length follows from 
equations (2.9) and (2.12) and reads 



l{k) = 1(2 



+ 5- 



(2.13) 



The same method allows a calculation of the mean-square deviation of the length 
of animals with area k. The generating function 



C{z)^f2c,z' ^. = E^(^-l)(Ll) 



(2.14) 



fc=i / 
is given by the second derivative of F{z, t) ai t = 1, 

d'^F{z,t) 



C{z) 



2zF'\z) + 2z^F^{z) 



(2.15) 



The coefficients of the series F^{z) are needed to obtain C^. They are calculated in 
the appendix with the following result: 



Thus we have 



fe=3 



F^z) = J: Ft%^ Ff = (^±iM - - . 



5 



2{k-l)F,-kF,_, 



and, using (2.12), 



5/c^ 1 



25 



k{k-l) 



k-2 



+ k\ 13-4 



+ 7 



(2.16) 



(2.17) 



(2.18) 



-r2 , 



The mean-square deviation A.P{k) = P{k) — I (k) follows from (2.13) and (2.18) 



2.3. Mean height 

The mean height of an animal with area k living on / stairs is taken as the ratio k/l, 
i.e. it is measured from the staircase. Thus the mean height of animals with area k is 
given by 



_ EZi kl-'F,, _ B, _^k(k-l 



(2.19) 
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The absorption/extraction identity can be used to write 



k-l + 1 



k-l + 1 /k-l 
I W-1 



k + 1 /k-l 
I W-1 



k-l 
l-l 



1>1, 



from which we deduce: 

kfk-l 
IW-1 



/c+ 1 



k-l + l\ (k-l 
I j^U-1 



I > 1 



Inserting (2.21) into the definition of i?^ in (2.19), we obtain 



Bu = 



k 



k + 1 
k 

k+\ 



i>i 



i-i 



(2.20) 



(2.21) 



(2.22) 



where the last equation follows from adding and substracting the term Z = in the 
first sum. According to the combinatorial definition of the Fibonacci numbers in (2.7) 
we have 



k 



k + 1 



iFk+2 + F,-1) 



(2.23) 



and finally 



h{k) 



k + 1 



(2.24) 



2.4- Asymptotic behaviour 

The generating function of the Fibonacci numbers in (2.8) can be written as the partial 
fraction expansion 



where ^ 1.61803 is the golden ratio and = 1 — 0. Thus the Fibonacci numbers are 
given by 

i^fe = ^(0'-0')- (2.26) 

Generally the number of animals with size k behaves at large size as Ck~^\^, 
where C is a constant amplitude, A is the inverse of the critical fugacity and 9 is 
the exponent governing the critical behaviour of the generating function G{z, 1) ~ 
{z^ — zy~^ when z — > z^_. From (2.7) and (2.26) we deduce the critical parameters: 



(2.27) 
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According to equations (2.13) and (2.26), the mean length behaves asymptotically as 

~m = ^ k + 0(1) = k + 0{l). (2.28) 

For the mean-square length we have 

Pik)='-±^e+'-^^k + Oil) (2.29) 

whereas 

W=(i±«!.^ + 5(^.^0(1). (2.30) 

It is easy to verify that the leading 0(/c^) contributions in (2.29) and (2.30) are the 
same. Thus the mean-square deviation is 0{k) and given by: 

k 



AP(/c) = ^ + 0(l). (2.31) 

The behaviour of the mean height follows from equations (2.24) and (2.26). It 
tends to a constant value: 

h{k) ^l + 4>^ + 0{k-^) = + 0{k-') . (2.32) 



3. Double staircase 

Next we consider column-convex animals on a double staircase as shown in figure 2. 
The connectivity rules for neighbouring columns are the same as above for a single 
staircase and we count the different stackings with at least one square on the lowest 
central stair. 

3.1. Number of animals 

The generating function for the number ^ of animals with area k living on I 
successive stairs, 

oo 

G{z,t)= J2 Gk^iz't' (3.1) 

k,l=l 

can be expressed using the generating function for animals in a single staircase in 
equation (2.3) as 

^2 

G{z,t) = zt + t- [1 + F{z,t)f. (3.2) 

1 — z 

The first term corresponds to a single square on the central stair. When there are two 
or more squares on the central stair (first factor in the second term) the animal may 
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Figure 2. Different stackings of five squares on a double staircase arranged by 
number of occupied stairs. 



either stop or continue to grow on one or two staircases (second factor in the second 
term). Using (2.2), the generating function may be rewritten as 



G{z, t) = zt + [l + F{z, t)] [F{z, t) - zt] = (1 - zt) F{z, t) + F'^{z, t) . 
Accordingly, we have 



G 



_ [zH']F{z,t) - [z''-H'-^]F{z,t) + [zH']F\z,t) 

k-l 
l-l 



.1-2 

Using the relation [13, p 169] 

— q\ fk — I — p + q 



k — l\ sr^ fp — q\ fk — p — I + q 

p,q ^ 



I 



1 



E 



^q — 1 J \ I — q — 1 
the double sum in (3.4) can be reduced to 



k-l + l 
l-l 



i-i 



EE 

9=1 P 

yielding 



p — q \ / k — p — I + q\ fk — I + 1 



G 



q-1 



k-l 
l-l 



l-q-1 



l-l 

E 

q=l 



l-l 



{l-l] 



k-l + 1 
l-l 



where the last expression follows from the addition/induction relation (2.5). 
The total number of animals with area k reads as 

«.^E«.^E<^^tV^E(L1). 



(3.3) 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



According to (2.9), the first term is equal to whereas the second follows 

from (2.7). Using the expression of Af, in (2.12) and the recursion on the Fibonacci 
numbers, we obtain 



G, 



A,+, - 2F,_, 



k{F, + 2F,_,) + 4F, - 5F,_, 



(3.9) 
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3.2. Mean length 

As above the mean length of animals with area k is defined as 

^'^'^itf^^ •^-i:<Vit^)-i:<L1). (3.10) 

The are the coefficients of the generating function 



dt 
fc=i 



= -zF(z) + il-z) A{z) + 2F(z)A(z) (^-^l) 

= -zF{z) + ^i^^ F\z) + 2 F\z) 
z z 

which follows from (3.3), taking into account the expression of A{z) given in (2.10). 
Thus, according to (3.11), (2.11) and (2.16), we have 

A, = [z']Aiz) = Ff_\- 2Ff )+ Fi%+ 2 [f^^- Ff 

^ 5fc^F,_, + /c(19F, - 7F,_,) + 6F, - 25F,_, ^^"^^^ 

25 

The mean length follows from equations (3.9), (3.10) and (3.12) 

T..^ _ 1 ^^'^k-i + - 7F,_,) + 6F, - 25F,_, 

^^"5 fc(i^. + 2F,_J+4F,-5F,_, " ^'"^'^ 

The mean-square deviation can be determined, proceeding as above for the single 
staircase. It involves a lengthy calculation of the generating function F'^{z), which 
can be obtained in the same way as F^{z) in the appendix. The asymptotic behaviour 
is the same as for the single staircase as shown in the discussion. 

3.3. Mean height 

The mean height, measured from the staircase, takes the form 



I 



I 



According to (2.7), the first sum is equal to F^_|_j^. Using (2.5) and (2.19), the second 
sum can be rewritten as 

^^^k/k-l\ ^k/k-l + l\ ^kfk-l\ k „ 

Et(,_2)=Ei( i_i )-Ei(,_ij=ifcTT^'"-^- 

Collecting these results and taking (2.23) into account, we obtain 

2k 2k 

= kF^+i + ^ (^fe+i + - 1) - ^q:^ (3F,+i + - 1) 

(3.16) 

= (.+ l)(. + 2) [^^^-^)^^^^ + ^^^ + ^)^^-^]- 
Finally, the mean height follows from (3.9) and (3.16) and reads 

^ _ 5fe k\F^ + Ffe_J + fe(F, - Ffe_,) + 6F,_, 

(/c+l)(/c + 2) MF, + 2F,_,) + 4F,-5F,_, ' ^ 
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3.4- Asymptotic behaviour 

According to (3.9) and (2.26), the asymptotic number of animals with area /c on a 
double staircase is such that 

9 = -l \ = (j). (3.18) 

From (3.13) we deduce that 

m = ^ + 0(1) = k + 0(1) (3.19) 



whereas (3.17) leads to 

2 + 



Xk) = 5 + 0{k-^) = + 0{k-^) . (3.20) 



4. Discussion 

The behaviour of the number of animals F^^ as a function of Z for A; » 1 can be 
obtained by expandind In ^ to second order near its maximum using the Stirling 
approximation. This leads to a Gaussian distribution: 



V27ffc 



Fk,l - exp 



(/ ^k)f 
2AP{k) 



(4.1) 



Here l(k) and AP(A;) are the leading contributions to the mean length and the mean- 
square deviation as given in (2.28) and (2.31), respectively. The prefactor follows from 
the normalization to ~ (p'^/^/d. 

The same behaviour is obtained in the double staircase for which, according 
to (3.7), G^i ~ IF^i — \k)Fj^i. Only the normalization differs since, according 
to (3.9), Gfc ~ k{Ff. + 2Ff._^)/h ~ k(j)^/b. Hence we have 



2AP{k) 



(4.2) 



Thus, to leading order, l{k) and AP{k) are the same for both problems. 
Furthermore the Gaussian distribution leads to h{k) — kl/l{k) ~ k/l{k) and we 
obtain the relation h{k)l{k) — k, valid to leading order too. These conclusions are in 
agreement with the exact results obtained in subsections 2.4 and 3.4. 

The animals are strongly anisotropic. Their length, measured along the staircase, 
scales as l{k) ~ /c'^n and their transverse size as h{k) ~ k'^ with z^y — vz = 1 and 
V — Thus the anisotropy exponent z is infinite. 

In this paper we have studied column-convex animals in staircases with steps of 
unit height. The problem can be generalized by considering staircases with steps 
of constant arbitrary height. Then Fibonacci numbers are replaced by generalized 
Fibonacci numbers. 
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Appendix. Calculation of the coefficients of F^{z) 

Starting from the partial fraction expansion for F{z) in (2.25) we hav: 



5^5 



(1-0^)3 (1-0^)3 {l-(j)z){l-^z) l-^z 



(A.l) 



Making use of the series expansion 



k + 2 



(1_^)3 Z_.V k 



E 



X 



(A.2) 



fc=0 



with X = (f)z and a; = 02; in the two first terms and, in the third, taking into account 
the relation 



{1 - (pz){l - 4>z) ^\-z-z^ 
which follows from the expressions of and in (2.25), we have 

1 



(A.3) 



F\z) 



Furthermore, since 



5\/5 



1 V^{7\ 



k=0 



(A.4) 



ik Ik r^fc 



= [z'] 



= V5K 



1-02; 1-02, 
we can rewrite (A.4) as 

3 _ (fc + l)(fc + 2) p . _ 3 



(A.5) 



k=l 



(A.6) 



Thus, the coefficients of the series are given by 



(3)_ (A; + l)(A- + 2) 
10 

(/c + l)(fe + 2) 
10 



F.-^\z' 



(A.7) 



^k r^k+1- 



Using (2.11), it is easy to verify that = ^2^^^ = and Fg^^^ = 1 as expected since 
the series expansion of F{z) starts with z. 
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